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1.  Background 


As  part  of  an  effort  to  model  single-walled  graphene  carbon  nanotubes,  a  transformation  is  first 
described  which  maps  atom  locations  originally  on  a  planar  sheet  to  atom  locations  on  a 
cylindrical  nanotube.  The  mapping  is  parametrized  by  the  two  components  describing  the  chiral 
vector,  and  thus  is  valid  for  generating  the  atomic  locations  of  arbitrary  chirality  nanotubes. 
Attention  is  then  turned  towards  quantifying  the  dimensions  of  the  unit  cell,  again  parametrized 
by  the  chiral  vector’s  two  components.  Taken  together,  the  mapping  plus  the  unit  cell 
dimensions  are  used  to  generate  unique  atomic  locations  of  arbitrarily  long  carbon  nanotubes. 

A  second  mapping  is  described  which  generalizes  the  commonly  used  Born  rule  to  higher  orders. 
The  Bom  rule  is  an  effective  tool  for  linking  bulk  material  deformations  to  atomic 
displacements,  particularly  for  homogeneous  deformations.  For  nonhomogeneous  deformations, 
such  as  those  involving  nonzero  curvature,  the  standard  Born  rule  can  inaccurately  link  the  two 
scales.  The  higher  order  mapping  described  herein  allows  nonhomogeneous  deformations  to  be 
mapped  down  to  the  atomic  scale  to  an  arbitrary  degree  of  precision. 

The  two  mappings  described  herein  are  anticipated  to  be  important  developments  leading  to 
accurate  deformation  modeling  of  graphene  carbon  nanotubes  of  arbitrary  chirality:  the  first 
mapping  will  allow  for  the  generation  of  the  initial  atomic  locations  of  the  undeformed  carbon 
nanotube,  and  the  second  mapping  will  allow  subsequent  deformations  of  the  nanotube  to  be 
accurately  modeled  at  the  atomic  level. 


2.  Introduction 


The  Bom  rule  is  an  important  component  in  the  development  of  continuum  models  of  atomic 
lattices.  It  serves  as  an  elementary  kinematic  assumption  that  establishes  a  connection  between 
atomic  motion  and  continuum  gradients,  and  its  origins  can  be  found  in  the  work  of  Cauchy 
(1828a,  1828b,  1829a,  1829b),  who  hypothesized  that  the  gradients  of  the  deformation  agree 
with  the  motion  of  molecules  in  a  deforming  solid.  Later,  in  his  first  book  on  dynamics  of 
crystal  lattices.  Bom  generalized  the  hypothesis  of  Cauchy  by  allowing  for  the  separate 
relaxation  of  internal  microstmcture  (Born,  1915).  This  text  was  later  translated  into  English 
(Bom  and  Huang,  1954).  Accounting  for  internal  relaxation  has  important  implications  for  many 
types  of  crystal  stmctures  and  crystal/continuum  mechanics,  e.g.,  in  elastic  properties  of  crystals 
in  the  absence  or  upon  removal  of  centrosymmetry  (Cousins,  1978a,  1978b,  1979,  1982;  Cousins 
and  Martin,  1978),  in  deformations  involving  phase  transformations  (Ericksen,  1984;  Zanzotto, 
1992,  1996),  in  contributions  to  the  elastic  properties  from  n-body  interactions  (Martin,  1975a, 
1975b,  1975c),  local  instabilities  in  the  deforming  Bravais  lattice  (Eriesecke  and  Thiel,  2002), 
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and  in  generalized  schemes  for  transitioning  atom  mechanics  to  continuum  mechanics  (Friesecke 
and  James,  2000),  among  others. 

In  general,  a  crystalline  lattice  L  in  its  reference  (unloaded)  configuration  is  composed  of  a 
number  of  interpenetrating  Bravais  lattices  whose  points  are  given  by 

X  =  M  'a,.  +p^ ,  with  i  =  1,  2,  3,  M‘  g  Z,  (1) 

where  X  are  the  lattice  points,  a,  are  the  linearly  independent  lattice  vectors  or  Bravais  base 
vectors,  and  are  the  shift  vectors  for  the  inner  atoms.  For  N  +  1  atoms  in  the  basis,  the  index 
k  mns  from  0  to  N.  The  base  vectors  a,  generate  a  so-called  skeletal  Bravais  lattice  S  czi 
where  each  point  in  S  along  with  the  k- shifted  atoms  compose  a  cluster,  or  the  so-called 
microstructural  motif  (see  Zanzotto,  1992).  When  the  crystal  is  undeformed,  the  base  vectors  are 
maximal  or,  namely,  they  represent  the  smallest  possible  periodic  unit  of  the  crystal.  When  the 
crystal  is  subjected  to  a  deformation,  the  base  vectors  a,  are  no  longer  maximal  in  general.  For 

this  reason,  fet  us  call  R  —\X  —  M  'a,  ,M '  e  z}  the  maximal  skeletal  lattice  of  L  .  Let  D  be 

the  set  of  all  allowable  deformation  gradient  tensors  F  and  let  R  ’  be  the  maximal  skeletal 
lattice  in  the  deformed  configuration  generated  by  the  set  of  lattice  vectors  a' ■ .  Then  the  formal 

definition  of  the  Bom  rule  is 

fR  =  R  ,  i.e.,  a'=m/Fa,  with  {ml)e  GL{3,Z),  for  all  FgD  .  (2) 

The  notation  GL(3,  Z)  comes  from  group  theory  referring  to  the  classical  linear  group  of 
invertible  3x3  matrices  with  integer  coefficients,  which  effectively  offers  the  ability  to  select 
larger  sublattices  that  are  maximal  in  their  representation  of  crystal  periodicity  in  the  context  of 
deformations. 

For  completeness,  it  is  noteworthy  to  mention  the  hypothesis  of  Cauchy,  which  is  more 
restrictive  than  Born’s  and  states  that  for  all  Fg  D  ,  [  ‘is  the  deformed  configuration 
corresponding  to  the  undeformed  [  ,  such  that  fL  =i  ‘.  The  so-called  Cauchy  mle,  or  Cauchy- 
Born  mle,  says  that  the  atomic  motion  agrees  with  the  gross  deformation  of  the  crystal,  whereas 
the  Bom  mle  says  that  only  the  motion  of  the  skeletal  lattice  agrees  with  the  gross  deformation 
and  the  microstmctural  motif  is  free  to  relax  into  its  appropriate  equilibrium  configuration.  It 
follows  then  that  when  the  base  vector  of  the  skeletal  lattice  corresponds  to  the  vector  connecting 
nearest  neighbors  in  a  single  species  crystal,  the  Cauchy  and  Born  hypotheses  are  equivalent. 

In  general  deforming  crystals,  the  Bom  mle  for  a  maximal  lattice  R  can  be  surprisingly 
restrictive  despite  the  freedom  to  choose  a  maximal  skeletal  lattice  in  L  ,  In  particular,  and  as  is 
discussed  further,  the  linear  nature  of  the  mle  prevents  it  from  being  applied  to  certain  types  of 
reduced  dimension  problems  undergoing  nonhomogeneous  deformations,  phenomena  involving 
diffusion,  lattice  shuffles,  or  local  instabilities  (Cousins  and  Martin,  1978;  Zanzotto,  1996; 
Friesecke  and  Thiel,  2002;  Arroyo  and  Belytschko,  2002).  Another  possible  approach  for 
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connecting  atomic  motion  and  continuum  deformations  is  described  by  Zanzotto  (1992)  through 
the  so-called  weak  Bom  rule,  which  allows  greater  flexibility  in  the  selection  of  the  sublattices 
but  is  still  subject  to  similar  requirements  as  the  Born  rule.* 

In  spite  of  some  of  its  limitations,  the  Bom  rule  continues  to  be  applied  in  advanced 
computational  methods  with  success  (see  Tadmor  et  ah,  1996;  Shilkrot  et  ah,  2002).  In  general, 
the  mle  has  been  found  to  yield  excellent  results  in  three-dimensional  (3-D)  bulk  or  two- 
dimensional  (2-D)  planar  materials.  The  concern  in  this  report,  however,  is  in  its  application  to 
reduced  dimensions  such  as  thin  films  embedded  in  3-D  space.  The  motivation  for  such 
problems  arises  naturally  in  contemporary  problems  of  atomistic  nanotubes  and  thin  films,  where 
the  nanotube  geometry  can  be  constmcted  from  a  planar  graphene  sheet.  The  key  issue  in  the 
modeling  of  such  applications  is  the  kinematic  hypothesis  inherent  in  the  Bom  mle,  which 
assumes  that  some  level  of  homogeneity  exists  in  the  deformation  map.  While  the  concept  of 
homogeneous  deformation  is  quite  standard  for  a  3-D  body,  it  is  not  the  same  for  generalized 
geometric  surfaces.  These  objections  were  noted  by  Arroyo  and  Belytschko  (2002)  who 
recognized  a  fundamental  distinction  between  (1)  tangents  of  the  undeformed  and  deformed 
bodies  that  are  mapped  from  one  to  the  other  via  the  deformation  gradient  F  and  (2)  the  chords 
created  by  lines  connecting  two  atoms  in  a  curved  manifold.  The  chords  are  of  critical 
importance  because  they  determine  the  relative  energy  of  an  atom  to  its  neighbor.  For 
generalized  geometric  surfaces,  chords  are  not  mapped  correctly  by  the  deformation  gradient. 

Arroyo  and  Belytschko  recognized  the  need  to  map  chords  correctly  for  generalized  geometric 
surfaces  and  therefore  proposed  an  alternative,  but  approximate,  method  to  handle  Bom  mle- like 
kinematics  on  curved  manifolds  using  the  so-called  exponential  Bom  mle.  At  its  core  is  an 
exponential  map  that  maps  the  tangent  of  a  manifold  back  into  the  manifold  itself,  thereby 
intrinsically  parameterizing  the  surface  with  the  local  tangents.  The  scheme  serves  as  a 
correction  step  following  the  application  of  the  traditional  Born  rule.  Although  effective,  it 
requires  knowledge  of  the  surface  via  geodesics  which,  in  general,  involves  the  solution  of  a 
nonlinear  system  of  ordinary  differential  equations  whose  unknowns  are  the  parametric 
coordinates  of  the  geodesic.  For  the  special  case  of  a  surface  that  forms  a  uniform  cylinder,  the 
calculation  is  straightforward  and  can  be  done  analytically  when  accompanied  by  a  geometric 
approximation  regarding  atom  placement  on  a  curved  manifold  (Arroyo  and  Belytschko,  2003). 
The  end  result  is  a  membrane  theory  based  only  on  the  geodesic  curvature  of  the  cylinder,  in  the 
form  of  a  correction  that  accurately  maps  one  generalized  cylindrical  surface  into  another.  It 
would  be  useful  for  nanotube  applications  to  explicitly  consider  the  nanotube  chiral  indices 
(Harris,  1999)  into  the  mapping  scheme. 


Our  definition  of  the  Born  rule  differs  slightly  from  that  of  Zanzotto  (1992).  We  presume  that  the  Born  rule  pertains  only  to 
the  mapping  of  the  skeletal  Bravais  lattice  and  not  to  the  selection  of  the  sublattice  itself.  The  “requirements”  mentioned  here  are 
in  reference  to  the  assumptions  of  homogeneous  deformations  in  the  absence  of  instabilities,  diffusion,  and  preservation  of  a 
Euclidean  metric. 
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In  this  report,  we  propose  an  alternative  method  to  the  Bom  mle  and  exponential  Bom  mle  for 
general  manifolds.  In  doing  so,  we  explicitly  consider  nanotube  indices  for  specific  application 
to  carbon  nanotubes  and  develop  the  closed-  form  expressions  for  mapping  an  arbitrary  graphene 
sheet  into  a  cylindrical  configuration.  The  larger  aim  is  to  determine  the  two-axis  bending 
rigidity  and  torsional  rigidity  of  a  graphene  carbon  nanotube  of  arbitrary  chirality.  As  a  first  step 
toward  realizing  this  aim,  attention  is  initially  focused  on  developing  a  mapping  for  graphene 
sheets  of  arbitrary  chirality  which  take  Bravais  indices  {n,  m)  and  locate  atomic  positions 
{x,  y,  z)  on  the  surface  of  a  cylindrical  carbon  nanotube.  In  order  to  ensure  a  mapping  which  is 
one-to-one,  bounds  on  (n,  m)  are  developed  which  define  the  domain  of  the  so-called  unit  cell. 

The  second  part  of  this  report  details  the  proposed  mapping,  termed  herein  the  higher  order  Born 
mle,  intended  to  replace  the  standard  Bom  mle  when  linking  bulk  material  deformations  to 
atomic  deformations.*  In  contrast  to  the  exponential  Born  rule,  the  higher  order  Bom  rule 
eliminates  the  need  for  a  correction  step  and  is  considerably  quicker  (in  a  computational  sense) 
and  easier  to  implement.  Depending  on  the  order  to  which  the  mle  is  taken,  it  can  be  used  to 
yield  a  scale- linking  rule  of  arbitrary  precision. 


3.  General  Chirality  Graphene  Sheet  to  Nanotube  Mapping 


3.1  Mapping 

In  this  section,  a  mapping  'T  :  {n,m)  g  (x,  y,  z)  g  is  developed,  taking  integers  n  and  m, 

indicating  the  usual  Bravais  indices  for  a  graphene  sheet,  to  position  (x,  y,  z)  on  a  manifold  N, 
representing  the  surface  of  a  graphene  carbon  nanotube.  Specifically,  consider  the  planar 
graphene  sheet  of  figure  1  in  which  aj  and  a2  denote  Bravais  base  vectors,  C  =  Na^  +  Ma^ 
denotes  the  chiral  vector,  T  =  Aa^  -i-  Ma2  denotes  the  translation  vector  normal  to  C ,  and  0 

indicates  the  angle  C  makes  with  the  reference  Cartesian  frame  (X,  Y,  Z)  with  unit  vectors  i 
and  j  .  The  length  of  the  chiral  vector  can  be  found  in  any  standard  reference  on  lattice 

geometry, 

C  =  ||C||  =  l^-^N^  +  NM  +  M^  ,  (3) 

o 

where  4  denotes  the  length  of  the  basis  vectors  aj  and  a2 ,  which  for  graphene  is  0.246  A  .  The 
length  of  T  determines  the  unit  cell,  defined  to  be  the  smallest  rectangle  defined  by  C  and  a 
translate  of  C  such  that  all  four  corners  of  the  unit  cell  coincide  with  an  atomic  lattice  point. 


We  should  also  note  that  at  the  time  of  writing,  another  paper  (Sunyk  and  Steinmann,  2003)  on  a  higher  order  gradient 
concept  appeared  that  discussed  the  application  of  the  Cauchy -Bom  mle  to  inhomogeneous  in-plane  (2-D)  deformations  of 
atomistic/continuum  problems.  The  work  presented  herein,  however,  differs  in  the  generalization  of  the  Bom  rule  for  arbitrary 
precision  and  its  application  for  reduced  dimensions  problems  in  a  full  3-D  sense. 
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Figure  1.  Lattice  geometry  for  the  graphene  sheet. 

This  length  is  also  well  documented,  see  for  example  Harris  (1999),  and  is  given  by, 

T  = 

where  denotes  the  highest  common  divisor  of  N  and  M  and  z  denotes  any  integer.  Any 
lattice  point  on  the  graphene  sheet  can  be  located  with  a  position  vector  r  and  lattice  indices  n 
and  m,  or  chiral  and  translational  coordinates  c  and  t, 

r  =  na^+ma2  =  ce^  +  te^,  (5) 

where  unit  vectors  along  C  and  T  are  denoted  by  and  e, ,  respectively. 

Based  on  this  description  of  the  graphene  sheet,  a  mapping  can  be  developed  which  describes  the 
positions  of  atoms  on  a  carbon  nanotube  formed  by  rolling-  up  the  graphene  sheet  along  C ,  such 
that  the  two  end-points  of  C  occupy  the  same  point  on  the  nanotube,  and  such  that  T  indicates 
the  tubular  direction.  The  development  of  this  mapping  starts  first  with  the  simple  map 
taking  points  (A,T,Z)g  to  (x,y,z)&R^, 

x  =  d>i(X,T,Z)  =  rsin(X/r),  y  =  d>2(X,T, Z)  =  T  ,  z  =  d>3(A,T,Z)  =  rcos(A /r)- r ,  (6) 

which,  when  applied  to  points  on  the  planar  manifold  fl,  IT  =  {(X,T,0)  1 0  <  X  <  r,0<Y  <  /}, 
yields  points  on  a  new  cylindrical  manifold  X  with  radius  r  and  axis  {(O,  y,  r)  1 0  <  y  <  /},  as 
shown  in  figure  2. 


II  II  I — ;; - ; — ; - 77  (  ^/3C/d„  N-M^Sidf, 

It\\  =  IJN^ +NM'+M'^  =\  ^  2  ",  " 

"  "  [Sc/{3d„)  N-M  =  3zdi, 


(4) 
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Figure  2.  Mapping  taking  planar  manifold  fl  to  cylindrical  manifold  X . 

The  mapping  just  described  is  that  required  to  roll-  up  a  graphene  sheet  along  chiral  vectors  given 
by  (N,  0).  To  roll-up  the  sheet  along  any  general  chiral  vector  (N,  M),  must  be  applied  first  to 
points  referenced  to  a  coordinate  system  aligned  with  the  chiral  and  translational  directions,  i.e., 
(c,  t,Z)  G  ,  thereby  rolling  up  an  angled  sheet  fig  (angle  to  the  horizontal  measured  by  0)  into 

an  angled  cylinder  Xg .  Following  this  initial  mapping,  a  second  mapping  dj''  (corresponding  to 
a  coordinate  transformation)  must  be  applied  to  transform  the  new  points  on  X  g  referenced  to 

the  chiral  and  translational  axes  to  points  referenced  to  the  (X,  Y,  Z)  system.  The  composition 
mapping  is  denoted  by  'F *  =  o  d> :  (c,  f,  Z)  g  (v,  y ,  z)  g  R^  and  is  given  as 

X  =  'F*(c,f,Z)  =  rcos(0)sin(c/r)-fsin(0), 

y  =  'F2*(c,f,Z)  =  rsin(0)sin(c/r)-l-fcos(0). 


z-'^^{c,t,Z)  =  rco^{clr)-r,  (7) 

where  d**^ :  (c,  t,  Z)  {x,  y,  z)  is  represented  in  matrix  form  by 


X 

COS0 

-sin  0 

O' 

c 

y 

•  = 

sin  0 

COS0 

0 

t 

» 

9 

(8) 

z_ 

0 

0 

1 

Z 
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and  where  a  superscript  R  indicates  this  mapping  holds  for  points  in  and  not  in  Z^,  as 
desired. 

The  desired  mapping  'T  follows  closely  from  'T*  .  Noting  that  the  Bravais  base  vectors  for 
graphene  are  given  as 


1  ^2 


• 

2  ^0^  2  ’ 


(9) 


expressions  for  coordinates  (c,  t)  can  be  developed  using  equations  5,  8,  and  9, 


c-Iq - cos(0)  +  7o - sin(0),  t  =  -4 - sin(0)+  4 - cos(0),  (10) 


while  expressions  for  0  and  r  can  be  developed  in  terms  of  chiral  indices  N  and  M, 


0  =  sin  ' 


1  VSM 

+  NM  +M  - 


J.  t-Q 

r  =  — - 


1 


+  A^M  +  M‘ 


V  J 

2  ,  .  /  „  ^  d3  • 


71 


(11) 


The  desired  map  ^  :{n,m\N ,M)&  Z  {x,y,z)  &  R  is  then  arrived  at  after  further 
manipulation 


X  ='Pj(7i,7w; 7V,M)  =  — 4.^4{7V^+Mf  +M^)-3M^  sin 
471 


n-\ — 
2 


7t 


|74{7V^  +  Mf  +M^)-3M^  +  ^twiM 


N^  +  NM  +M^ 


+ 


4  3(277  +  777 )M "  - 3mM^4(N^  +NM  +  M^)-3M‘ 
8  N^  +  NM  +  M^ 


='1^2(77, 777;  A^,  M)  = — 4 


471 


sin 


^  777  ^ 

77  H - 

2 


71 


|J4(a^"+A^M+M")-3M"  +  ^TimM 


N^+NM+M^ 


4  (277  +  m)SM  - mSM^4(N^  +  NM  +M^)-3M‘ 
8  N^  +  NM  +  M^ 


-.J4(n^  +  NM+M^)-3M'^ 


^  ^  777^ 


■=^^(n,m-,N,M)=—ljN^  +  NM+M 


2n 


cos 


71 


77  H - 

2 


^4(a^"+A^M  +  M")-3M"  +^mnM 


N^+NM  +  M^ 


V 

-^^N^+NM  +M'^ 
2k 


(12) 
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where  arguments  after  the  ;  indicate  chiral  parameters.  An  example  of  the  map  is  shown  in 
figure  3  for  a  (9,  5)  chiral  vector. 


Figure  3.  Cylinder  (blue)  generated  using  mapping 

'F  :(n,  m;  A?, M)g  (x, y,z)e where 

{N  =  9,  M  =5)  chirality  has  been  chosen  and 
where  the  range  of  indices  used  to  generate  the 
figure  are  —  5  <  n  <5 ,  0<m  <8-  The  red  plane 
indicates  the  locus  of  (X,  Y,  Z)  locations  using 
the  given  indices  and  before  applying  the  map. 


3.2  Unit  Cell 

Application  of  'P  maps  points  from  a  plane  to  a  cylinder,  but  without  proper  choice  of  the 
indices  ranges  for  (n,  m),  results  in  an  overdetermination  of  the  cylindrical  manifold — i.e.,  the 
map  is  not  one-to-one,  as  more  than  one  point  on  the  plane  is  mapped  to  the  same  point  on  the 
cylinder.  To  prevent  this  occurrence,  bounds  on  the  indices  are  developed  next  which  ensure  a 
domain  on  the  plane  which  is  mapped  one-to-one  onto  the  cylinder. 

Attention  is  focused  on  the  unit  cell  of  figure  1  in  order  to  develop  the  indices  ranges.  For  chiral 
vectors  with  associated  0  <  60°,  the  minimum  index  for  which  a  lattice  point  lies  in  the 
unit  cell  is  0,  while  the  minimum  index  is  determined  by  N' .  Expressions  for  N'  and  M' 
can  be  developed  using  the  translation  vector  and  its  magnitude, 

T  =  r(- sinG  i-i-cos6  j)=  Aki -b  ,  (13) 


with  equation  9.  The  desired  expressions  are 


The  maximum  index  is  simply  N,  while  the  maximum  index  is  determined  by  the  m 

index  of  a  vector  formed  by  the  addition  of  T  and  C — i.e.,  the  top- right  comer  of  the  unit  cell.  In 
order  to  determine  this  index,  it  is  convenient  to  introduce  the  dual  basis  vectors  and  such 
that  a'  •  a^  =  5 where  the  center  dot  denotes  the  Cartesian  inner  product, 


Then, 


1  1 

r. 

2  1 

a  =  — 

1 — j 

,  a  — 

— j 

lo 

3 

(  ) 

3 

(  y 

"^max  =(T+C)-a" 


(15) 


IT 


COS0  -I-  M  . 


(16) 


Note  that  use  of  equations  4  and  11  in  equations  14  and  16  yields  indices  bounds 
(^min  ’  ^max  ’  "^min  ’  ^max )  which  UTC  Only  functions  of  the  chiral  vector’s  N  and  M. 


With  bounds  identified  for  the  indices  m  and  n,  two  conditions  must  still  be  met  on  any  pair  of 
indices  (n-  <n<  ,m„.  <  m  <  )  to  ensure  their  existence  inside  of  the  unit  cell: 

0  <  c  <  C  and  0<t<T  where  it  is  recalled  that  c  and  t  are  chiral  and  translational  coordinates 
corresponding  to  (n,  m)  and  are  given  by  equation  10  while  C  and  T  are  the  chiral  and  translation 
vector  magnitudes  given  by  equations  3  and  4.  Given  a  chiral  vector  {N,  M),  a  single  nested  loop 
can  be  easily  programmed  using  the  minimum  and  maximum  bounds  on  n  and  m  plus  the 
conditional  checks  to  generate  a  list  of  index  pairs  of  atoms  located  in  the  unit  cell.  The  mapping 
'T  given  in  equation  12  can  then  be  applied  to  locate  the  atoms  on  the  nanotube  manifold  N. 
Results  of  this  procedure  are  shown  in  figure  4  where  unit  cells  and  nanotubes  for  several 
common  chiralities  are  shown. 


4.  On  a  Higher  Order  Born  Rule 


The  Bom  mle  is  commonly  employed  to  determine  atomic  lattice  positions  from  the  bulk 
material  deformation  gradient  F.  Specifically,  letting  R,  denote  the  position  vector  of  the  atom 
in  the  continuum  before  deformation,  and  letting  r,  denote  the  position  of  the  same  atom  after 
deformation;  the  Bom  mle  relates  the  two, 

r,=F-R,,  (17) 

where  F  referenced  to  Cartesian  coordinates  (Xi,  X2,  X^)  in  the  undeformed  configuration  and 
ix\,X2,  X3)  in  the  deformed  configuration  can  be  stated  as 
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Figure  4.  Results  of  procedure  for  determining  the  indices  in  the  unit  cell  (left  column)  and  subsequent  mapping 
(right  column)  for  chiral  vectors  (a)  N  =  5,  M  =  5,  (h)  N  =  6,  M  =  3,  and  (c)  N  =  9,  M  =  0. 
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(18) 


F 


dr 


®e. 


e,.  0e 


j  ■ 


As  discussed  by  Arroyo  and  Belytschko  (2002),  the  deformation  gradient  maps  points  on  the 
tangent  space  associated  with  the  undeformed  configuration  to  points  on  the  tangent  space 
associated  with  the  deformed  configuration,  and  therefore  inadequately  maps  deformations  with 
nonzero  curvature.  They  address  this  deficiency  by  applying  the  Bom  mle  as  stated  in 
equation  17,  and  then  correcting  ri  through  the  exponential  mapping.  Note  that  this  correction 
follows  the  application  of  the  Born  mle.  Furthermore,  by  their  own  account,  the  implementation 
of  the  exponential  map  is  “not  straight  forward”  and  requires  additional  simplifications.  Herein, 
a  higher  order  Bom  rule  is  proposed  which  improves  upon  the  standard  Bom  rule  by  improving 
the  initial  mapping  from  R,  to  r„  obviating  the  need  for  a  follow-up  correction,  in  contrast  to  the 
approach  of  (Arroyo  and  Belytschko,  2002). 

As  opposed  to  equation  17,  the  differential  element  dR  maps  to  Ar  exactly  through  the  mle 

dr=FdR,  (19) 


which  can  be  integrated  once  to  obtain, 

r  =  j  F  ■  dR .  (20) 

Further  integration  of  the  right-hand  side  requires  an  expression  for  F  =  F(R) .  It  is  chosen  to 
expand  F  in  a  Taylor  series  about  R  =  0  without  loss  of  generality. 


3F 

F(R)  =  F(0)+-^ 

oR 


R  +  6>(R") 


R=0 


(21) 


Substitution  of  equation  21  into  20,  and  denoting  the  constant  third-order  tensor 
yields  the  desired  higher  order  Born  mle, 

r=F(0)  R  +  K  Jr  dR  +  6>(R") . 


dF 

dR 


R=0 


as  K  , 


(22) 


Comparison  of  equation  22  with  the  standard  Born  mle  (equation  17)  reveals  the  underlying 
assumption  of  the  Born  rule — that  the  deformation  gradient  is  constant  over  all  R  in  which  it  is 
applied.  This  is  an  adequate  approximation  for  homogeneous  or  small-  gradient  deformations 
like  expansion,  contraction,  or  shear,  but  may  be  inadequate  for  high-  gradient  deformations  or 
for  bending- like  deformations  with  nonzero  curvature. 

Next,  an  example  deformation  mapping  is  exercised  in  order  to  demonstrate  the  utility  of  the 
higher  order  Born  rule.  Consider  again  the  mapping  illustrated  in  figure  2,  taking  points 
(Aj,  Aj,  A3)e  to  (x^,  X2,x.^)e  , 
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=  $j(Xj,X2,  X3)  =  rsin(Xj  /r),  ^2  ^3)  =  ^2’ 

^3  =  ^>3(Xj,X2,X3)  =  rcos(Xj  I  r)-r . 

The  deformation  gradient  for  this  map  is 

cos(Zi/r)  0  O' 

0  1  0 
-sin(Zi/r)  0  0 


(23) 


(24) 


while  the  only  nonzero  component  of  K  is  K3jj  =  -1  / r .  Note  that  this  indicates  negative 
curvature  in  the  z-v  plane  with  radius  of  curvature  r  for  all  (Xj,  X2,  X3) ,  as  expected.  Note 
further  that  the  higher  order  Born  rule  (equation  22)  for  this  example  evaluates  as 

1  V  2 

r  =  Xiei+X2e2 - ^e3+6>(R"),  (25) 

r  2 


where  the  nonzero  63  coefficient  demonstrates  the  ability  of  the  higher  order  Bom  mle  to  more 

accurately  map  points  from  the  plane  to  the  cylinder.  If  desired,  further  accuracy  can  be  obtained 
by  taking  more  terms  in  the  expansion  (equation  21).  In  contrast,  the  standard  Bom  mle  for  this 
mapping  is  inadequate  and  would  only  map  the  plane  onto  itself. 


5.  Conclusions 


Two  important  mappings  have  been  developed  in  an  effort  to  accomplish  the  larger  goal  of 
determining  the  two -axis  bending  and  torsional  rigidity  of  graphene  carbon  nanotubes.  The  first 
mapping  allows  atomic  lattice  coordinates  to  be  generated  for  graphene  carbon  nanotubes  of 
arbitrary  chirality.  In  addition  to  the  mapping,  bounds  on  the  Bravais  indices  have  been 
developed  (thereby  describing  a  unit  cell)  which  ensure  the  mapping  is  one-to-one.  The  second 
mapping  developed  allows  subsequent  deformations  of  the  nanotube  to  be  linked  to  atomic 
displacements  using  the  newly  termed  higher  order  Bom  mle.  Although  the  new  mle’s  use  has 
been  demonstrated  in  the  context  of  graphene  carbon  nanotubes,  it  is  expected  to  be  equally 
valuable  in  the  study  of  mechanics  associated  with  any  generalized  geometric  surface  of  reduced 
dimension,  which  is  embedded  in  3-D  space. 
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